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The phase diagram at finite temperature and density is investigated in the framework of the 
Polyakov linear sigma model (PLSM) with three light quark flavors in the mean field approximation. 
It is found that in the PLSM, the three phase transitions, i.e, the chiral restoration of u, d quarks, 
the chiral restoration of s quark, and the deconfinement phase transition are independent. There 
exists two-flavor quarkyonic phase at low density and three-fiavor quarkyonic phase at high density. 
The critical end point (CEP) which separating the crossover from the first-order line in the PLSM 
model is located a.t{TE, i^e) = (188 MeV, 139.5 MeV). In the transition region the thermodynamic 
O . properties and bulk viscosity over entropy density ratio C,/ s are also discussed in the PLSM. 
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r-| . I. INTRODUCTION 

^ • It is widely believed that at sufficiently high temperatures and densities there is a quantum chromodynamics (QCD) 
D \ phase transition between normal nuclear matter and quark-gluon plasma (QGP), where quarks and gluons are no 

^ . longer confined in hadronsQ. The theoretical and experimental investigation of QGP is one of the most challenging 
problems in high energy physics. The main target of heavy ion experiments at the Relativistic Heavy Ion collider 

. (RHIC), the forthcoming Large Hadron Collider (LHC) and FAIR at GSI is to create such form of matter and study 
^ ■ its properties. 

, There are two different phase transition which are associated with two opposite quark mass limits. For Nf massless 
CN . quark flavors, the QCD Lagrangian possesses a chiral U{Nf)R x U{Nf)L = SU{Nf)H x SU{Nf)L x f/(l)y x 1/(1) a 
[ symmetry, here V = R + L, while A = R ~ L. However this symmetry does not appear in the low energy particle 
. . spectrum, it is spontaneously broken to the diagonal SU{Nf)v group of vector transformation by a non-vanishing 
I expectation value of the quark-antiquark condensate, (^rIl) 7^ 0. This process involves iV| Goldstone bosons which 
dominate the low-energy dynamics of the theory. The U(l)v symmetry is always respected and thus plays no role 
in the symmetry breaking pattern considered in the following discussion. The axial U{\)a symmetry is broken to 
Z{Nf)A by instanton effects [|||^. Consequently, one of the N'j Goldstone bosons becomes massive, leaving Nj — 1 
Goldstone bosons. The SU{Nf)]i x SU^Nf)^ x U{1)a group is also expficitly broken by the effects of nonzero quark 
masses. The N'j — 1 low energy degrees of freedom then become pseudo-Goldstone bosons. For M < Nf degenerate 
] quark flavors, an SU{M)v symmetry is preserved[QlP. It is generally believed that at sufficiently high temperature 
and density there should be a transition from ordinary hadronic matter to a chirally symmetric plasma of quark and 
gluons. At temperature of about 170 MeV, Lattice QCD calculations indicate that this symmetry is restored]^. The 
order of the phase transition seems to depend on the mass of the non-strange u and d quarks, and the mass of the 
strange quark rris, and at the temperature on the order of 200 MeV, heavier quark flavors do not play an essential 
role. 

On the contrary, in the heavy quark limit, QCD reduces to a pure SU{Nc) gauge theory which is invariant under 
a global Z{Nc) center symmetry, where the Polyakov loop serves as an order parameter|0, H, ^ |l^. In general, the 
Polyakov loop is a complex scalar field and it is related to the free energy of a static quark in the gluon medium, it 
vanishes in the confining phase means that the quark has infinite free energy and takes a finite value in the deconfined 
phase. It is linked to the Z{Nc) center symmetry of the SU{Nc) gauge theory, thus the confining phase is center 
symmetric, while the center symmetry is spontaneously broken in the deconfined phase. In the present of dynamical 
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quarks the center symmetry is explicitly broken. No order parameter is established for the deconfinement transition 
in this case, but the Polyakov loop still serves as an indicator of a rapid crossover towards deconfinement . On 
the other hand, the chiral transition has a well-defined order parameter in the chiral limit of massless quarks: the 
quark-antiquark condensate, ^Jip/qi) 7^ 0, which and its dynamical generation, is the basic element of the linear sigma 
model and the Nambu-Jona-Lasinio model. 

Recently, the effective chiral models for QCD with U{Nf)ii x U{Nf)L symmetry, such as the Nambu-Jona- 
Lasinio model(NJL) or, via bosonization, the linear sigma model (LSM) could be combined with the Polyakov 
loop which allows to investigate both, the chiral and the deconfinement phase transition. These models named as 
the Polyakov-Nambu-Jona-Lasinio model (the PNJL) JT2, O, 04), and the Polyakov linear sigma model (PLSM) 



or the Polyakov quark meson model (PQM)|T| ^ \l7l , are widely studied in recent years. In the PNJL and 
PLSM model, where quarks develop quasiparticle masses by propagating in the chiral condensate, while they cou- 
ple at the same time to a homogeneous background (temporal) gauge field representing the Polyakov loop dy- 
namics. The full QCD thermodynamics at zero and finite quark chemical potential have been investigated in 
Refs.|T|, 011, [T|, ^, |2l|, H 11, H, H 1^, 1^, H, 1^, l^. However, most studies focus on the PNJL model 



with two and three quark flavors. As we know both of the NJL model and the LSM model for the phenomenology 
of QCD can be parameterized to describe equally well the vacuum structure at T = = McV, but these models 
treat the contribution of the Dirac sea differently, in the NJL model it is included explicitly up to a momentum 
cutoff A, while in the sigma model this contribution is renormalized out. So that the extended PLSM model which is 
augmented with the Polyakov loop has the benefit of renormalizability. In this paper, calculations for the thermody- 
namic potential and the phase diagram will be performed at the mean-field level as that of the PNJL model, however, 
eventually, such model will be investigated beyond mean-field by the loop expansion |3l) , especially by using the C JT 
effective potential p^. 

In present paper, we extend the Polyakov linear sigma model with two quark flavors JTsI , [l^ to three quark flavor [pT| 
to investigate the phase diagram at finite chemical potential. Moreover, since the study of the equation of state in 
QCD with the 2-1-1 quark flavors by the lattice QCD simulations with almost physical quark masses have recently been 
reported in Ref [^ at zero chemical potential but finite temperature, it is interesting to describe QCD thermodynamics 
in the PLSM with three quark flavors, and compare with the recent results of the lattice simulations at zero chemical 
potential, especially, the essentials of QCD thermodynamics around the critical temperature Tc- The outline of the 
paper is as follows: in the next section wc introduce the PLSM model with three quark flavors, some symmetry 
breaking patterns in the vacuum are briefly discussed, and parameters arc fixed. In Sect. Ill, after obtaining the 
effective potential in the mean field approximation, we explore the phase diagram at T — /i plane. Section IV is devoted 
to derive the thermodynamic properties of the system at zero chemical potential, i.e. the pressure, the equation of 
state, the square of the speed of sound, the specific heat, the trace anomaly of the energy-momentum tensor and the 
bulk viscosity, all these thermodynamical observables are compared with the lattice QCD data. At the end, we give 
discussions and summary in Sec.V. 



II. THE MODEL 

Following Ref. |^ , we introduce a generalized Lagrangian of the linear sigma model for Nf = 3 quarks and A^c = 3 
color degrees with quarks coupled to a spatially constant temporal background gauge field representing Polyakov loop 
dynamics (the Polyakov-linear-sigma model or the PLSM in short), the Lagrangian reads [ p^ 

C = Cchiral -l^i4>,4'* ,T) (1) 

where we have separated the contribution of chiral degrees of freedom and the Polyakov loop. 

The chiral part of the Lagrangian, Cchirai — + ^m^^^, of the SU{3)b. x SU{3)l symmetric linear sigma model 
with three quark flavors consists of the fermionic part 

= II V^/(*7^^M - 9Ta{<ya + il57^a))i>f (2) 
/ 

and the purely mcsonic contribution 

£,„ = Tr(a^$^a^$-m2$t$)_ Ai[Tr($t$)]2 
-A2Tr($t$)2 + c[Det($) + Dct($t)] 

+Tr[i?($ + $^)], (3) 

the sum is over the three flavors (f=l,2,3 for u, d, s). In the above equation we have introduced a flavor-blind Yukawa 
coupling g of the quarks to the mesons and the coupling of the quarks to a background gauge field A^^ — S^^qAq via 
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the covariant derivative = 9^ — *j4^j. The $ is a complex 3x3 matrix and is defined in terms of the scalar ctq and 
pseudoscalar tTq meson nonets, 

•P^Taiaa + ina). (4) 
The 3x3 matrix H breaks the symmetry explicitly and is chosen as 

H = Taha, (5) 

where ha are nine external fields. The Ta ~ Aa/2 are the generators of the U{3) symmetry, are the Gell-Mann 



matrices with Aq = y ^1. The Ta are normalized to TT{TaTi,) ~ 5ab/'2 and obey the U{3) algebra with [Ta, Th] ~ ifabcTc 

and {Tq, Tfc} = datcTc respectively, here fabc and date for a,b,c= 1, 8 are the standard antisymmetric and symmetric 
structure constants of SU{^) group and 

faba = 0, dabo = \[^^°-'^' 

In Eq.(H), is the tree-level mass square of the fields in the absence of symmetry breaking, Ai and A2 are 
two possible quartic coupling constants, and c is the cubic coupling constant. The terms in the first line of Eq.(|^) 
are invariant under [/(3)i^ x C/(3)l = [/(3)v x U{i)A transformations, the determinant terms are invariant under 
5t/(3)ij X SU{3)l — SU{3)v x SU{3)a, but break the U{1)a symmetry explicitly, these terms arise from the U{1)a 
anomaly of the QCD vacuum. The last term in Eq.(^) breaks the axial and the SU{3)v vector symmetries explicitly. 
A non- vanishing vacuum expectation value for $, ($) = Ta&a, breaks the chiral symmetry spontaneously. Because the 
parity is not broken in the vacuum, there are no non- vanishing vacuum expectation values for fields tt^. The patterns 
of explicit symmetry breaking have been described in detail in Ref. ^ for the ?7(3)y x [/(3)a model, in this work 
we will constrain our study on the case of explicit chiral symmetry breaking with U{\)a anomaly. Since ($) must 
carry the quantum numbers of the vacuum, only the fields CTq corresponding to the diagonal generators of C/(3) can 
be nonzero. The same holds for ha, then the diagonal components /iq, h^ and ft-s of the explicit symmetry breaking 
term could be nonzero. Furthermore due to the fact that the masses of the up and down quarks are approximately 
equal, rriu — rud, and the strange quark mass is larger than niu and nid, in the following discussion, we restrict 
our study to /iq 7^ 0, = and hs ^ 0. In this case the SU{3)v x SU{3)a symmetry is explicitly broken to SU{2)v 
isospin symmetry, and the U{1)a subgroup of the U(S)a symmetry is already explicitly broken by the instanton. 
Traditionally, the linear sigma model does not have quark degrees of freedom, the nonzero quark masses correspond 
to a term of the form Eq.(^) in the lagrangian, where the matrix H is proportional to the quark mass matrix, the 
fields H can be determined from the vacuum values for the pion and kaon masses, as well as the pion and kaon decay 
constants, the values for the constants of proportionality could be fixed at bare quark masses for u, d and s quarks, 
for example ruq = 10 MeV, rus = 150 MeVjSSt. This means that the parameters in the purely mesonic Lagrangian 
carries the information of the bare quark mass, in order to avoid the double counting, such model do not introduce 
bare quark masses in Lagrangian (^. 

The quantity U (0, (f>* , T) is the Polyakov-loop effective potential expressed by the dynamics of the traced Polyakov 
loop 

= {Tt,L)/N„ 4>* = (Tr.Lt)/^^. (7) 
The Polyakov loop L is a matrix in color space and explicitly given by 



L{x) — Vcxp 



drAilx, t) 







(8) 



with P = 1/T being the inverse of temperature and A4 = iA'^. In the so-called Polyakov gauge, the Polyakov-loop 
matrix can be given as a diagonal representation Jl2| . The coupling between Polyakov loop and quarks is uniquely 
determined by the covariant derivative -D^i in the PLSM Lagrangian in Eq.(|l|), and in the chiral limit, this Lagrangian 
is invariant under the chiral flavor group, just like the original QCD Lagrangian. The trace of the Polyakov-loop, (f> 
and its conjugate 0* can be treated as classical field variables in this work. 

The temperature dependent effective potential L({(f>, 4>* ■, T) is used to reproduced the thermodynamical behavior of 
the Polyakov loop for the pure gauge case in accordance with lattice QCD data, and it has the Z{?,) center symmetry 
like the pure gauge QCD Lagrangian. In the absence of quarks, we have (p ~ cf)* and the Polyakov loop is taken 
as an order parameter for deconfincment. For low temperatures, U has a single minimum at = 0, while at high 
temperatures it develops a second one which turns into the absolute minimum above a critical temperature Tq, and 
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the Z(3) center symmetry is spontaneously broken. In this paper, we will use the potential h({(f>, 4>*,T) proposed in 
Rcf . |13|| , which has a polynomial expansion in (j> and cj)* : 



T4 2 6 



with 

62(T) = ao + ai(^^j+a2(^^j + ^3 [Y ) ' ^^^^ 

A precision fit of the constants a^, bi is performed to reproduce the lattice data for pure gauge theory thermodynamics 
and the behavior of the Polyakov loop as a function of temperature. The corresponding parameters are 

ao = 6.75, ai = -1.95, 02 = 2.625, 

ag = -7.44, 63 = 0.75, 64 = 7.5. (11) 

The critical temperature To for deconfinement in the pure gauge sector is fixed at 270 MeV, in agreement with the 
lattice results. 

The chiral part of the model involves seven parameter m^, Ai, A2, c, ho, and two unknown condensates 
(To and (Tg, which we take from p^ , where they were fixed to reproduce some physical quantities in the hadronic 
sector. The vacuum condensates cto and ag £^re members of the scalar {Jp = 0^) nonet and both contain the up-down 
quark condensate and the strange quark condensate. It is more convenient to convert the condensates into a pure 
up-down quark condensate and strange condensate, this is achieved by an orthogonal basis transformation from the 
original octet-singlet basis {ao, ag) to the up-down quark condensate ax and the strange quark condensate ay basis 
(for simplicity, we ignore the bar over the up-down quark condensate ax and the strange quark condensate ax ), 

^5-0 + -^^8, (12a) 

1 



-V - ^cro - Y gCTg. (12b) 

As a result, when chiral symmetry breaks spontaneously, the fields ($) acquires a non-vanishing vacuum expectation 
value, the constituent quark masses for u, d and s are defined to be 

rriq = gax/2 for u,d, (13) 
TTis — gay 1^/2 for s, (14) 

and the light quark sector decouples from the strange quark sector. At zero temperature and zero chemical potential, 
the parameters of the Lagrangian are fixed in a way that these masses agree with the observed value of pion mass 
TOtt = 138 MeV and the most commonly accepted value for sigma mass mo- = 600 MeV. The values of the condensates 
are determined from the pion and kaon decay constants by means of the partially conserved axial-vector current 
relation (PCAS), 

f^x = /tt, CTy = ;^(2/a' - /tt), (15) 

and the decay constants, /tt = 92.4 MeV and = 113 MeV. The Yukawa coupling is usually fixed by requirement 
that the constituent quark mass in the vacuum, for u and d quarks, is about 1/3 of the nucleon mass, which gives 
g ~ 6.5. The strange constituent quark mass is predicted to be rUg ~ 433 MeV. The six parameters to^, Ai, A2, c, 
hx and hy are set for to„ = 600 MeV with the axial U{1)a anomaly: = (342.52MeV)^ , Ai = 1.40, A2 = 46.48, 
c = 4807.84 MeV, hx = (120.73MeV)3 and hy = (336.41MeV)3||. 



III. MEAN FIELD APPROXIMATION 



The standard approach for dealing with the thermodynamics of variable particles is via the grand canonical ensemble. 
Let us consider a spatially uniform system in thermodynamical equilibrium at temperature T and quark chemical 
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potential ^f(f = u,d,s). In general, the grand partition function reads 

Z = Trexp[-(7T;- ^ ^^fNf)/T] 



f—u,d^s 

" " " (16) 



f—u^d,s 

where j^=i j^^'^ dt Jy d^x and V is the volume of the system. The u and d quark chemical potentials is 

/"« = + Y + -MY, = 2" 3^^' " "3~ ^ 3^^ ^ ^ 

where fiB, fJ-i and fiy are the baryon, isospin and hypercharge chemical potentials respectively. In general, the three 
quark chemical potentials are independent, but in the following discussion, we assume the SU{2)v isospin symmetry 
by neglecting the slight mass difference between an u~ and d~ quark. Then the light quark chemical potentials 
become equal, /i^ = /i„ = ^cj, and there is no effect of the isospin density since the isospin chemical potential /i/ is 
zero. Also for simplicity and confronting our results with the other models, such as LSM and PNJL models with three 
quark flavors, we will only consider symmetric quark matter and define a uniform chemical potential fi = fiq = /i, for 
Hy = 0. 

We evaluate the partition function in the mean-field approximation similar to p4[ | |3^ . Thus we replace the meson 
fields by their expectation values in the action. In other words, we neglect both quantum and thermal fluctuations of 
the meson fields. The quarks and antiquarks are retained as quantum fields. The integration over the fermions yields 
a determinant which can be calculated by standard methods[0. This generates an effective potential for the mesons. 
Finally, we obtain the thermodynamical potential density as 

n{T,^lf) = = Uia,,ay) +Ui(t>,(f>*,T) + n^^, (18) 



= -2TiV, / ^{ln[l + 3(0 + re-(^'-^)/^)xe-(^«-^)/^ + e-3(B,-M)/T] 



with the quarks and antiquarks contribution 

d^p 

+ln[l + 3(0* + 0e-(^'+^)/^) X e-(^'+^)/^ + e-3(iJ,+M)/T]| 
-2TNs I 7?|r{ln[l + 3(0 + 0*6-'^-^)/^) x e'^^-'^)/^ + e-3(^^-M)/T] 

J y/''^) 

+ln[l + 3(0* + 0e-(^=+^)/^) x e-(^=+^)/^ + e-3(^^+'')/^]}. (19) 



Here, Nq = 2, Ng = 1, and Eg = ^/jp + is the valence quark and antiquark energy for u and d quarks, for strange 



quark s, it is E's = -^/p^ + mf, and niq, rris is the constituent quark mass for u, d and s in Eqs.(|l3|)([l4|). The purely 
mesonic potential is 



U{cr^,ay) = ^{crl+cr'fj) - K<t^ - hyffy - '^'^I'^y + y^x^^y + ^(2Ai + Aa)^;^ + ^(Ai + A2)cr^. (20) 



Minimizing the thermodynamical potential in Eq.(|l8|) with respective to ct^., ay, and 0*, we obtain a set of equations 
of motion 

5f2 _ ^ - - - (21) 

da^ ' day ' 90 ' 90* 

The set of equations can be solved for the fields as functions of temperature T and chemical potential /i, and the 
solutions of these coupled equations determine the behavior of the chiral order parameter ax, ay and the Polyakov 
loop expectation values 0, 0* as a function of T and fi. 



IV. T- fj, PHASE DIAGRAM 



We now explore the phase diagram of the Polyakov linear sigma model at finite temperature and density. There 
are three different critical temperatures, i.c, for the chiral phase transition of u, d quarks, for the chiral phase 
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transition of s quark, and T^f for the deconfincment phase transition. In order to locate the critical temperatures, we 
define the temperature derivatives of the condensates: cr^ = da^/dT, a'y = dcjy/dT, </>' = d(j}/dT and cj)*' = dcj)* /dT. 
In the following, we will show the double peak structure of the temperature derivatives of the condensates and explain 
explicitly how we determine the critical temperatures. 



A. Phase transitions at zero density 

We firstly investigate the phase transitions at zero density. In Fig.0(a), the temperature dependence of the chiral 
condensates Ux for u, d quarks, the chiral condensate ay for s quark and the Polyakov loop expectation value 0, 0* 
at ^ = MeV is shown in relative units. Here and in the following discussion, the chiral condensates arc normalized 
by their zero-temperature value: axo = 92.4 MeV and Cyo = 94.5 MeV for u, d and s quarks, respectively. The 
temperature derivatives of the condensates a'x,(y'y,<i>' and </>*' are shown in Fig||(b). 

The temperature behavior of the chiral condensates and Polyakov loop condensate shows that the system experiences 
a crossover at zero chemical potential. The temperature derivative of the chiral condensate a^, for u and d quarks has 
only one peak at T ~ 200 MeV. The temperature derivative of the chiral condensates CTj, for s quark has two peaks, 
the left peak coincides with the peak of at T ~ 200 MeV, the right peak shows up at T ~ 260 MeV. For vanishing 
chemical potential, we find = 0*, and ~ 1.1 for T — > oo. The temperature derivatives 0' and 0*' shows one peak 
and one pseudo-peak, and the peak coincides with the peak of cr^ at T ~ 200 MeV. 

Now we show how we determine the critical temperature of the crossover. When there is only one peak in the 
temperature derivative of the condensate, the location of the peak gives the critical temperature T^. When there 
are two peaks in the temperature derivative of the condensate, for chiral phase transition, the critical temperature 
is determined by the peak temperature corresponding to a ^ ,y{T) / a ^ ,y{T = 0) < 1/2, and for deconfinement phase 
transition, the critical temperature is given by the peak or pscudo peak temperature corresponding to 0(T)/0(T — > 
c») > 1/2. This definition of the critical temperature is different from that defined in Ref . [p^ . 

The temperature derivatives of the condensates in Fig.|l| (b) show that the deconfinement phase transition happens 
at a higher critical temperature ~ 220 MeV, and the chiral restoration occurs at a smaller critical temperature 
~ 203.5 MeV. This result is similar to that obtained in the two-flavor Polyakov linear sigma model Ref.||l^, and 
different from that in the three-flavor Nambu-Jona-Lasinio model with the Polyakov loop jlj], where they find the 
simultaneous crossovers around ~ 200 MeV. The critical temperature for the chiral restoration of the strange quark 
is at Tf ~ 260 McV. 
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FIG. 1: (a) The normalized chiral condensate (Ta;, Oy and the Polyakov loop 0, as a function of temperature for = 
MeV. (b) Temperature derivatives of the normalized chiral condensate Ua;, dy and the Polyakov loop 0, 0* as a function of 
temperature at /i = MeV. The Polyakov variable and dy are scaled by a factor of 6. 
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B. Phase transitions at finite density 



For fixing the chemical potential at 200 MeV, the order parameters and their temperature derivatives as functions 
of temperature are shown in Fig.H (a) and (b), respectively. 

It is found that the temperature derivative of the chiral condensate for u and d quarks still has only one peak, 
at the peak temperature ~ 173 MeV of cr^, the chiral condensate jumps from 0.8 to 0.15, which indicates that the 
chiral phase transition for the u and d quarks is of first order phase transition. 

The temperature derivative of chiral condensates Oy for s quark has two peaks, the left peak coincides with the 
peak of < at ~ 173 MeV, the right peak shows up at Tf ~ 240 MeV. 

The values of the Polyakov loop <^ and (/>* are different for nonvanishing chemical potential, due to the fact that the 
free energies of quark and antiquarks are different in the finite chemical potential When [i increases, the left 

peak moves down to a smaller temperature, which is always coincident with T^, and the right pseudo-peak at zero 
chemical potential develops into a real peak and stays almost at the same temperature ~ 220 MeV as that for 
/Lt = 0. 



Here are several remarks on the phase transitions. The left peak of a'^ and 
of (T^, which reflects the effective coupling among a^^ and 



* ) always coincides with the peak 
However, the coupling among the condensates 
does not necessarily induce the simultaneous phase transitions. From the above analysis, we can see that in general 
cases, the three phase transitions, i.e, the chiral restoration of u, d quarks, the chiral restoration of s quark, and the 
deconfinement phase transition are independent in the Polyakov linear sigma model. However the lattice result shows 
the simultaneous phase transitions of chiral restoration and deconfinement. it should be very interesting to investigate 
how to realize the simultaneous phase transitions in the Polyakov linear sigma model, and it deserves further efforts 
on understanding the correlation between the chiral phase transitions and deconfinement phase transition. 
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FIG. 2: (a) The normalized chiral condensate a^^ ay and the Polyakov loop (j), (f>* as a function of temperature for /i = 200 
MeV. (b) Temperature derivatives of the normalized chiral condensate (Jx, <Jy and the Polyakov loop 0, (p* as a function of 
temperature at fi — 200 MeV. The Polyakov variable and ay are scaled by a factor of 6. 



C. T — ji phase structure and the location of the critical end point 



We evaluate the chiral phase transitions of u, d and s quarks and deconfinement phase transition at finite temperature 
and finite density, and show the T — fi phase structure of the Polyakov linear sigma model in Fig.H (a). 

From Fig.0 and Fig.^, we see that the chiral phase transition for s quark and the deconfinement phase transition 
for the Polyakov loop are always crossover, for two light flavors, there is a crossover in the low density region and a 
first-order phase transition in the high-density region, and in the middle exists a critical end point (CEP). In order 
to locate the critical end point, we define the quark number susceptibility Xq — d^fl/d^^, which should be divergent 
at the critical end point. In Fig. ^ (b), we plot the quark number susceptibility as a function of the temperature 
for different chemical potentials. For the Polyakov linear sigma model, the result shows that the critical end point 
is around (Te,^£;) = (188 MeV, 139.5 MeV), which is close to the lattice result {Te^he) = (162 ± 2 MeV, he = 
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120 ± 13 McV) |40j. For the linear sigma model without the Polyakov loop, the critical end point is located at 
{Te,He) - (92.5MeV, 216 MeV). The critical chemical potential he in PLSM is much lower than that in the PNJL 
model with three quark flavors where the predicted critical end point is ^e > 300 MeVp^, |2^ . 

The chiral phase transition for the strange quark and the deconfinement phase transition in the T — ^ plane 
are shown in Fig.|| (a) by the dash-dotted line and dotted line, respectively. It is found that with the increase of 
chemical potential, the critical temperature for strange quark to restore chiral symmetry decreases. However, for 
the deconfinement phase transition, with the increase of chemical potential, the deconfinement critical temperature 
keeps almost a constant around 220 MeV. It can be seen that in the Polyakov linear sigma model, there exists two- 
flavor quarkyonic phase pl[ at low density, where the m, d quarks restore chiral symmetry but still in confinement, and 
three-flavor quarkyonic phase at high density, where the w, d, s quarks restore chiral symmetry but still in confinement. 
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FIG. 3: (a) The T — jj, phase diagram in the Polyakov linear sigma model, (b) the quark number susceptibility Xg ioi different 
chemical potentials as function of T. 



V. THERMODYNAMIC PROPERTIES OF THE POLYAKOV LINEAR SIGMA MODEL 

In order to investigate the influence of the Polyakov loop on the equilibrium thermodynamics, we calculate several 
thermodynamic quantities. All information of the system is contained in the grand canonical potential which is given 
by in Eq.(|l8|), evaluated at the mean field level. The entropy density is determined by taking the derivative of 
effective potential with respect to the temperature, i.e, 

s = ^dn{(f>)/dT. (22) 

As the standard treatment in lattice calculation, we introduce the normalized pressure density p which is normalized 
to vanish at T = /i = and the energy density e as 

p^-n, £ = -p + Ts. (23) 

The equation of state p{e) is an important input into hydrodynamics. The square of the speed of sound is related 
to p/e and has the form of 

2 ^ dp ^ S ^ ^ 

' de Tds/dT a' ^ ' 

where 

a = de/dT, (25) 



is the specific heat. At the critical temperature, the entropy density as well as the energy density change most quickly 
with temperature, thus one expect that should have a minimum at T^. The trace anomaly of the energy-momentum 
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tensor T'"' 



A 



2^4 



3p 



y4 



(26) 



is a dimensionless quantity, which is also called the " interaction measure" . 

In Fig.^^, we show the trace anomaly, energy density, pressure density, equation of state parameter, specific heat 
and sound velocity square as functions of the scaled temperature T jT^. 

The result of the trace anomaly of the energy-momentum tensor T^^ at zero density is shown in Fig^(a) in 
comparison with lattice data for A^,- = 6 . It is shown that the trace anomaly approaches the conformal value at 
high temperature in both the linear sigma model and the Polyakov linear sigma model, which agrees with the lattice 
result. The trace anomaly shows a peak around in the linear sigma model, and the height is only half of that 
of lattice result. After including the Polyakov loop, the peak of the trace anomaly appears at a higher temperature 
1.5 Tc comparing with the lattice result, where the peak shows up at around 1.1 Tc- We will show later that the 
appearance of the peak in the trace anomaly is not related to the phase transition, but the change rate of the trace 
anomaly is related to the phase transition. 

The result of the trace anomaly of the energy-momentum tensor at different chemical potentials as function of the 
scaled temperature T/Tc is shown in Fig.^(b). It is observed that the trace anomaly changes smoothly around 
in the case of zero chemical potential when the system experiences a crossover, but changes sharply in the case of 
fist order phase transition at high density. From Figj^, we can see clearly that the behavior of the trace anomaly at 
Tc resembles that of the energy density and the entropy density at Tc. We will show latter that the change rate of 
the energy density and/or entropy density at Tc determine the critical behavior of other thermal quantities like the 
specific heat, the sound velocity square. 




T- 



■ • n=0 MeV 
■ n=80 MeV 
-H=139.5 MeV 
-H=16Q MeV 




( b ) 

FIG. 4: (a) The scaled trace anomaly (e - 3p)/r'* for = MeV. The Polyakov linear sigma model prediction (solid line) and 
the linear sigma model prediction (dashed line) are compared with N f = 2 + 1 lattice QCD data for Nr = 6. Lattice data taken 
from Ref.[B3|. (b) The scaled trace anomaly (e — Zp)/T'^ in the Polyakov linear sigma model for different chemical potentials 
as functions of T /Tc. 



Fig.^ shows the pressure density over energy density p/er, which is represented in terms of equation-of-state (EOS) 
parameter, at zero density and finite density, respectively. We observe that the pressure density over energy density 
increases with temperature and saturates at high temperature. Both the linear sigma model and the Polyakov linear 
sigma model give very similar results at high temperature, the pressure density over energy density p/ e saturates at 
a value smaller than 1/3. Another common feature of the p/e in the linear sigma model and the Polyakov linear 
sigma model is that there is a bump appearing at low temperature region, which is also observed in the lattice result. 
Around the critical temperature Tc, the pressure density over energy density p/e shows a downward cusp. However, 
the minimum value of the p/e around Tc is 0.2 in the linear sigma model, which is much larger than the result from 
the Polyakov linear sigma model and the lattice QCD data. For the Polyakov linear sigma model, the minimum of 
p/e around Tc is 0.075, which is consistent with the lattice QCD data |3^. When the chemical potential increases, 
from Fig.^ (b), we can see that the minimum of the p/e around Tc decreases. 

We show the specific heat as a function of the scaled temperature T/Tc in Fig.|^. It it shown that there is a sharp 
peak arising at Tc both in the Polyakov linear sigma model and linear sigma model. At high chemical potential when 
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FIG. 5: (a) The scaled pressure density p/T* in the Polyakov linear sigma model for different chemical potentials as functions 
of T/Tc- (b) The scaled energy density e/T'* in the Polyakov linear sigma model for different chemical potentials as functions 
of T/Tc 
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FIG. 6: (a) The equation-of-state parameter w{T) — p{T)/e{T) for fi — MeV. The Polyakov linear sigma model prediction 
(solid line) and the linear sigma model prediction (dash line) are compared with Nf =2+1 lattice QCD data for A^r = 6. 
Lattice data taken from Ref.|^^. (b) The equation-of-state parameter w{T) = p{T) / e{T) in the PLSM for different chemical 
potentials as functions of T /Tc- 



the phase transition is of first order, the specific heat diverges. From the definition of the specific heat = de/dT, 
it is easy to understand that the appearance of the peak is due the fast change of the energy density at the critical 
temperature T^. The other issue presented in the temperature behavior of in Fig.^is that there is a small second 
peak at higher temperature, which is more visible in the Polyakov linear sigma model than that in the linear sigma 
model. This result is consistent with the PLSM model with two quark flavors [p^ . 

Fig.^ shows the sound velocity square as a function of the scaled temperature. In conformal field theories 
including free field theory, the sound velocity square is always 1/3. In both the linear sigma model and the Polyakov 
linear sigma model, the saturates at a value smaller than 1/3 at high temperature. However, near the critical 
temperature Tc, a downward cusp shows up. At zero density, the minimum of the sound velocity square at Tc is 
around 0.05 and 0.16 in the Polyakov linear sigma model and the linear sigma model, respectively, and the value of 
Cl at Tc in the Polyakov linear sigma model is close the lattice result. When the chemical potential increases, the 
minimum of at Tc decreases, and it approaches in the case of first order phase transition. 
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FIG. 7: (a) The specific heat C„ as a function of the temperature for fi = MeV. The solid line denotes the Polyakov linear 
sigma model prediction and the dashed line denotes the linear sigma model prediction. (b)The specific heat Cv in the PLSM 
for different chemical potentials as functions of T/Tc. 
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FIG. 8; (a) The sound velocity square as a function of the temperature for ^ = MeV. The solid line denotes the Polyakov 
linear sigma model prediction and the dashed line denotes the linear sigma model prediction, (b) The sound velocity square 
Cs in the PLSM for different chemical potentials as functions of T/Tc. 

VI. THE BULK VISCOSITY OVER ENTROPY DENSITY RATIO IN THE POLYAKOV LINEAR 

SIGMA MODEL 



The bulk viscosity is related to the correlation function of the trace of the energy- momentum tensor 0^' : 

C = \ hm - r dt f d^re-* ([^^^(x), (0)]) . 



(27) 



According to the result derived from low energy theorem, in the low frequency region, the bulk viscosity takes the 
form of [Hi El 
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9wo I dT T4 
1 



9 Wo 



{-16e + 9TS + TCy + 16|e„|} 



(28) 
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with the negative vacuum energy density Sy = ily = Q{(f>)\T=o, and the parameter loq = i^a{T) is a scale at which the 
perturbation theory becomes vahd. From the above formula, we can see that the bulk viscosity is proportional to the 
specific heat Cy near phase transition, thus C/s behaves as l/C^ near Tc in this approximation. 
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FIG. 9: (a) The bulk viscosity over entropy density ratio (/s as a function of the temperature for fi — MeV. The solid line 
denotes the Polyakov linear sigma model prediction and the dashed line denotes the linear sigma model prediction. Lattice 
data taken from Ref. ||]. (b) The bulk viscosity over entropy density ratio ^/s in the PLSM for different chemical potentials 
as functions oiT/Tc. 



In Fig.^ (a) and (b), we plot the bulk viscosity over entropy density ratio C/s as a function of the temperature 
for zero chemical potential and finite chemical potential, respectively. From Fig.^ (a), it is shown that, at zero 
chemical potential /i = 0, the bulk viscosity over entropy density C^/s decreases monotonically with the increase of the 
temperature in both the Polyakov linear sigma model and linear sigma model, and at high temperature, C,/ s reaches 
its conformal value 0. In ||43[| , the bulk viscosity over entropy density of the three flavor system is extracted from 
lattice result, which is shown in Fig.^ (a) by the square. It is observed that C,/s in PLSM near phase transition is 
in very good agreement with the lattice result in P^ , i.c, it rises sharply near phase transition. From Fig.^ (b). 
which shows C,/ s as function of the scaled temperature T/Tc for different chemical potentials with = 0, 80, 139.5, 160 
MeV, We can see that when the chemical potential increases up to /i = 80 MeV, there is an upward cusp appearing 
in C/s right at the critical temperature T^- With the increase of the chemical potential, the upward cusp becomes 
sharper, and the height of the cusp increases. At the critical end point and when fi > fj.E for the first order phase 
transition, C/s becomes divergent at the critical temperature. 

The critical behavior of (^/s is determined by the shape of the trace anomaly or the energy density around the 
critical temperature Tc. In the case of first order phase transition, the energy density has a sudden change at the 
critical temperature, thus C/s diverges at Tc. In the case of crossover, C/s can exhibit different behaviors near phase 
transition: 1) When the energy density or entropy density changes slowly around Tc, ^/s monotonically but slowly 
rises up with the decrease of temperature. 2) When the energy density or the entropy density changes quickly near 
Tc, C/s monotonically but quickly rises up when the temperature decreases. 3) When the energy density or entropy 
density changes very quickly near Tc, an upward cusp of (^/s appears at Tc, i.e, C/s firstly rises quickly when the 
temperature decreases to Tc, then jumps fast and eventually rises up when the temperature drops further away from 
Tc. These behaviors have been analyzed in Ref. [Q for a toy model. 

It is expected that one can distinguish whether the system experiences a first order phase transition or a crossover 
from observables which are sensitive to the bulk viscosity at RHIC experiments. As mentioned in Refs.|^, ^ that a 
sharp rise of bulk viscosity near phase transition induces an instability in the hydrodynamic flow of the plasma, and 
this mode will blow up and tear the system into droplets. It would be interesting to investigate in more detail how 
different behaviors of (/s affect the observables. 

However, it is noticed that the results of bulk viscosity in this paper are based on Eq. (p8|), where the ansatz for 
the spectral function 
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has been used in the small frequency, and uiq is a scale at which the perturbation theory becomes valid. In our 
calculation, = 10 T, its magnitude at is in agreement with that obtained in ChPT for massive pion gas system 
in Ref . |^ . Qualitatively, the bulk viscosity corresponds to nonconformality, thus it is reasonable to observe a sharp 
rising of bulk viscosity near phase transition. Ref. ||4^ ] has investigated the correlation between the bulk viscosity 
and conformal breaking, and supports the results in Ref. p^, The sharp rising of bulk viscosity has also been 
observed by another lattice result and in the linear sigma model [Q. However, till now, no full calculation has 
been done for the bulk viscosity. The frequency dependence of the spectral density has been analyzed in Refs. ]5C| ] 
and and the limitation of the ansatz Eq.(p9|) has been discussed. From Eq. (psj), we sec that the bulk viscosity 
is dominated by at Tc- If Cy diverges at Tc, the bulk viscosity should also be divergent at the critical point and 
behave as t~°'. However, the detailed analysis in the Ising model in Ref. shows a very different divergent behavior 
( ~ with z ~ 3 the dynamic critical exponent and i/ ~ 0.630 the critical exponent in the Ising system. More 

careful calculation on the bulk viscosity is needed in the future. 



VII. SUMMARY AND DISCUSSION 



In this paper, we have extended the linear sigma model with three quark flavors to include certain aspects of 
gluon dynamics via the Polyakov loop. The PLSM model encodes two basic features that govern low energy QCD, 
spontaneous chiral symmetry breaking and confinement, and is a framework proposed to correctly interpret results 
from QCD thermodynamics and extrapolate to regions not accessible by lattice computations. 

Within the mean field approximation, we have studied the (T, /i) phase diagram of the Polyakov linear sigma model 
with the polynomial form of the Polyakov loop potential. It is found that in the linear sigma model with Polyakov 
loop, the three phase transitions, i.e, the chiral restoration of u,d quarks, the chiral restoration of s quark, and the 
deconfinement phase transition are independent and happen at different critical temperatures. It is found that in 
the PLSM, at low density, there exists two-flavor quarkyonic phase, where the u, d quarks restore chiral symmetry 
but still in confinement, and at high density there exists three- flavor quarkyonic phase where m, d quarks restore 
chiral symmetry but still in confinement. The linear sigma model with and without the Polyakov loop has a general 
feature: there is a crossover in the low density region and turn a first-order phase transition in the high-density region 
accompanied by a critical end point (CEP). Here we observe large discrepancies for the critical end point between the 
PLSM and LSM models. For the LSM model, the critical point is found at (Tb,//^) ~ (92.5, 216) MeV, while that of 
the PLSM model is at {Te,he) ^ (188, 139.5) MeV, which is close to the lattice result {TE,fiE) = {162±2 MeV, fiE = 
120 ± 13 MeV) [Q. The critical chemical potential fiE in PLSM is much lower than that in the PNJL model with 
three quark flavors where the predicted critical end point is at /x^ > 300 MeV[|2^, |2^ . 

In order to compare our results with the lattice QCD simulations and other models, at zero chemical potential but 
finite temperature, we have investigated the thermodynamic properties and bulk viscosity in the PLSM and LSM 
model. It is found that the inclusion of the Polyakov loop is necessary in order to quantitatively fit the lattice QCD 
result in the pure gauge sector. Our results in the PLSM and LSM models show that at critical temperature Tc, the 
trace anomaly A, the specific heat Cy show upward cusp at Tc. The ratio of pressure density over energy density 
p/e and the square of the sound velocity show downward cusp at Tc. These cusp behaviors at phase transition 
resemble lattice QCD results. We find that the PLSM model can reproduce all the thermodynamic and transport 
properties of the hot quark-gluon system near the critical temperature, p/e at Tc is close to the lattice QCD results 
0.075, the trace anomaly A ~ 6.6 at 1.3Tc is close to the lattice result of the peak value 7.7, the bulk viscosity to 
entropy density ratio (/s at Tc is around 0.2, which is also agreement with the lattice QCD result in Ref. p3|, psj . 

In this work the thermodynamics of the PLSM model have been compared to three-flavor lattice data, and it has 
been shown that the results of the model indeed agree with the lattice data, but it should be pointed that the used 
quark masses are different from each other, due to the fact that the lattice data with large quark masses while our 
model calculations physical values have been used. Furthermore because the model has several adjustable parameters, 
the numerical results of the critical temperature and the critical chemical potential of the phase transitions are, of 
course, parameter dependent. As above discussions for small temperatures and finite chemical potentials the chiral 
phase transition is probably of first-order while a crossover is expected at high temperatures and small chemical 
potentials, this suggests the existence of at least one critical cndpoint (CEP) where the first-order transition line in 
the phase diagram terminates. For the lattice simulations, although there are much progress have been achieved in 
the lattice studies for the QCD thermodynamics, at finite chemical potentials the fermion sigma problem is still a 
considerable obstacle, some lattice groups differ in their predictions. For the "CEP" problem, the location and even 
the possible existence in the phase diagram is still an open question ||5^[^4[, in such a manner, a comparison to the 
lattice data in Ref.Q is less meaningful. 

It is interesting to mention some further developments for this work, such as extensions beyond mean field theory 
by using loop expansion pi] or concerning both Polyakov loop and mesonic fluctuations |Q, and consideration 
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about matter at high baryon density for asymmetric quark matter, and such discussions are more relevant to the 
real experiment (RHIC and LHC) or observation (neutron stars), all of above extensions leave much room for further 
detailed investigations. Especially, at high chemical potential with asymmetric quark matter 7^ /id 7^ /£sj since in 
that case we can discuss the pion condensate and kaon condensate in the Polyakov linear sigma model |28[ |5^, ^t\ . 
Eventually, work in this direction is in progress. 
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